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Abstract. Gamma-ray bursts are thought to be the outcome of a cataclysmic 
event leading to a relativistically expanding fireball, in which particles are accel- 
erated at shocks and produce nonthermal radiation. We discuss the theoretical 
| predictions of the fireball shock model and its general agreement with observa- 

Q\ ■ tions. Some of the recent work deals with the collimation of the outflow and 

its implications for the energetics, the production of prompt bright flashes at 
wavelenghts much longer than gamma-rays, the time structure of the afterglow, 
its dependence on the central engine or progenitor system behavior, and the role 
of the environment on the evolution of the afterglow. 
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£ : I INTRODUCTION 

Gamma-ray bursts (GRB) have been studied in gamma-rays for over 25 
years, but except for rare and fleeting X-ray detections, until a few years 
ago there existed no longer-lasting detections at softer wavelengths. However 
in early 1997 the Italian-Dutch satellite Beppo-SAX suceeded in providing 
accurate X-ray locations and images that allowed their follow-up with large 
ground-based optical and radio telescopes. The current interpretation of the 
gamma-ray and longer wavelength radiation is that the progenitor trigger 
produces an expanding relativistic fireball which can undergo both internal 
shocks leading to gamma-rays, and (as it decelerates on the external medium) 
an external blast wave and a reverse shock producing a broad-band spectrum 
lasting much longer. 



A strong confirmation of the generic fireball shock model came from the 
correct prediction [43], in advance of the observations, of the quantitative 
nature of afterglows at longer wavelengths, in substantial agreement with the 
subsequent data [89,85,91,73,96]. The measured 7-ray fluences imply a total 
energy of order 10 54 (fi 7 /47r) ergs, where Af2 7 is the solid angle into which 
the gamma-rays are beamed. Collimation may indeed be present, evidence 
having been recently reported for this [34,18,12]. In any case, such energies 
are possible [44] in the context of compact mergers involving neutron star- 
neutron star (NS-NS ) or black hole-neutron star (BH-NS ) binaries, or in 
hypernova/collapsar models involving a massive stellar progenitor [56,68]. In 
both cases, one is led to rely on MHD extraction of the spin energy of a 
disrupted torus and/or a central BH to power a relativistic outflow. 

II THE GENERIC FIREBALL SHOCK SCENARIO 

Whatever the GRB trigger is, the ultimate result must unavoidably be an 
e ± , r y fireball, which is initially optically thick. The initial dimensions must 
be of order r min <ct var ~ 10 7 cm, since variability timescales are t var <10~ 3 s. 
Most of the spectral energy is observed above 0.5 MeV, hence the 77 — >• 
mean free path is very short. Many bursts show spectra extending above 1 
GeV, indicating the presence of a mechanism which avoids degrading these 
via photon-photon interactions to energies below the threshold m e c 2 = 0.511 
MeV. The inference is that the flow must be expanding with a very high 
Lorentz factor T, since then the relative angle at which the photons collide 
is less than T" 1 and the threshold for the pair production is diminished [25]. 
However, the observed 7-ray spectrum is generally a broken power law, i.e., 
highly nonthermal. In addition, the expansion would lead to a conversion of 
internal into kinetic energy, so even after the fireball becomes optically thin, it 
would be radiatively inefficient, most of the energy being kinetic, rather than 
in photons. 

The simplest way to achieve high efficiency and a nonthermal spectrum is 
by reconverting the kinetic energy of the flow into random energy via shocks 
after the flow has become optically thin [70]. Two different types of shocks 
may arise in this scenario. In the first case (a) the expanding fireball runs 
into an external medium (the ISM, or a pre-ejected stellar wind [70,39,30,78]. 
The second possibility (b) is that [71,54], even before external shocks occur, 
internal shocks develop in the relativistic wind itself, faster portions of the 
flow catching up with the slower portions. This is a generic model, which is 
independent of the specific nature of the progenitor. 

External shocks will occur in an impulsive outflow of total energy E Q 
in an external medium of average particle density n at a radius r<i ec ~ 
(3E /47rn m p c 2 ri 2 ) 1 / 3 ~ lO^-E^n" 1 / 3 ?^ 273 cm , and on a timescale tdec ~ 
r dec /(cT 2 ) ~3x lO 2 ^ 3 n" 1 ^" 8 / 3 s, where 77 = Y = 10 2 7? 2 is the final bulk 



Lorentz factor of the ejecta. Variability on timescales shorter than t dec may 
occur on the cooling timescale or on the dynamic timescale for inhomogeneities 
in the external medium, but generally this is not ideal for reproducing highly 
variable profiles [79] (see however [15]). However, it can reproduce bursts with 
several peaks [58] and may therefore be applicable to the class of long, smooth 
bursts. 

In a wind outflow [55], one assumes that a lab-frame luminosity L and 
mass outflow M Q are injected at r ~ r/ and continuously maintained over a 
time t w ; here 77 = L /M c 2 . In such wind model, internal shocks will occur 
at a radius [71] r dis ~ ct var rf ~ 3 x 10 14 t 1)ar ?7| cm, on a timescale t w ^> 
t var ~ rdis/ {erf) s, where shells of different energies Arj ~ rj initially separated 
by ct v (where t v < t w is the timescale of typical variations in the energy at 
ri) catch up with each other. In order for internal shocks to occur above 
the wind photosphere r p h ~ Max/(4:7rm p cT 2 ) = 1.2 x 10 14 L 5 3^ 3 cm, but 
also at radii greater than the saturation radius (so the bulk of the energy 
does not appear in the photospheric quasi-thermal component) one needs to 
have 7.5 x lO 1 ^ 5 ^ 1 / 5 <r]3 x lO 2 ^ 4 ^ 1 / 4 . This type of models have the 
advantage [71] that they allow an arbitrarily complicated light curve, the 
shortest variation timescale t var >,10~ 3 s being limited only by the dynamic 
timescale at r ; , where the energy input may be expected to vary chaotically. 
Such internal shocks have been shown explicitly to reproduce (and be required 
by) some of the more complicated light curves [79,33,61]. 

A potentially valuable diagnostic tool for the central engine of GRB is the 
power density spectrum (PDS). An analysis of BATSE light curves [3] indi- 
cates that the logarithmic slope of the PDS between 10~ 2 and 2 Hz is ap- 
proximately -5/3, and there is a cutoff of the average PDS above 2 Hz. Using 



FIGURE 1. Average power 
density spectrum Pf of sim- 
ulated bursts from internal 
shocks, compared with the ob- 
served PDS (thick line), using a 
square-sine modulated Lorentz 
factor and a cosmological dis- 
tribution satisfying the observed 
logN-logP (Spada, Panaitescu & 
Meszaros 1999). 

a simple kinematical model for the ejection and collision of relativistic shells 
[63,83] have calculated the light curves and PDS expected for a range of total 
burst energies and for a total mass ejected and bulk Lorentz factor distri- 
bution compatible with the internal shock scenario (Figure 1). The redshift 




distribution also affects the PDS, and the observed logN-logP relation is used 
as a constraint. For optically thin winds, a slope approaching -5/3 requires 
a non-random Lorentz factor distribution, e.g. with an asymmetrical time 
modulation so as to produce a larger number of collisions at low frequencies 
(see also [4]). A cutoff at high frequencies (~ 2 Hz) can be understood in 
terms of shocks which increasingly occur below the scattering photosphere of 
the outflow, or a deficit of energy in short pulses due to the modulation of the 
Lorentz factors favoring shocks arising further out. 

A significant fraction of bursts appear to have low energy spectral slopes 
steeper than 1/3 in energy [69,14]. This has motivated consideration of a 
thermal or nonthermal [35,36] comptonization mechanism, which can be put 
in the astrophysical context [23] of internal shocks leading to self-regulated 
pair formation. There is also evidence that the clustering of the break energy 
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FIGURE 2. Luminosity per decade xL x vs. x = hv/m e c 2 for two values of r\ = L/Mc 2 
and marginal (left) or large (rigt) pair compactness. T: thermal photosphere, PHC: pho- 
tospheric comptonized component; S: shock synchrotron; C: shock pair dominated comp- 
tonized component (Meszaros & Rees, 1999b). 

of GRB spectra in the 50-500 keV range may not be due to observational selec- 
tion [69,10,16]. Models using Compton attenuation [11] require reprocessing 
by an external medium whose column density adjusts itself to a few g cm~ 2 . 
More recently a preferred break has been attributed to a blackbody peak at 
the comoving pair recombination temperature in the fireball photosphere [17]. 
steep low energy spectral slope being due to the Rayleigh- Jeans part of the 
photosphere. In order for such photospheres to occur at the pair recombina- 
tion temperature in the accelerating regime requires an extremely low baryon 
load. For very large baryon loads, a related explanation has been invoked 
[87], considering scattering of photospheric photons off MHD turbulence in 
the coasting portion of the outflow, which upscatters the adiabatically cooled 
photons up to the observed break energy. These ideas have been synthesized 
[50] to produce a generic scenario (see Figure 2) in which the presence of 



a photospheric component as well as shocks subject to pair breakdown can 
produce steep low energy spectra and preferred breaks. 

Ill SIMPLE STANDARD AFTERGLOW MODEL 

The dynamics of GRB and their afterglows can be understood independently 
of any uncertainties about the progenitor systems, using a generalization of the 
method used to model supernova remnants. The simplest hypothesis is that 
the afterglow is due to a relativistic expanding blast wave, which decelerates 
as time goes on [43]. The complex time structure of some bursts suggests that 
the central trigger may continue for up to 100 seconds, the 7-rays possibly 
being due to internal shocks. However, at much later times all memory of the 
initial time structure would be lost: essentially all that matters is how much 
energy and momentum has been injected; the injection can be regarded as 
instantaneous in the context of the afterglow. As pointed in [70], the external 
shock bolometric luminosity builds up as L oc t 2 and decays as L oc t~( 1+<J ). 
Beyond the deceleration radius the bulk Lorentz factor decreases as a power 
law in radius, r oc r~ 9 oc r s/(1+29) , r oc t l ^ 1+2g \ with g = (3,3/2) for the 
radiative or adiabatic regime (in which pr 3 T ~ constant or pr 3 T 2 ~ constant). 

The synchrotron peak frequency in the observer frame is v m oc TB' / -f 2 , and 
both the comoving field B' and electron Lorentz factor 7 are expected to be 
proportional to T [39]. As T decreases, so will u m , and the radiation will move 
to longer wavelengths. For the forward blast wave, [53,31] discussed the pos- 
sibility of detecting at late times a radio or optical afterglow of the GRB. A 
more detailed treatment of the fireball dynamics indicates that approximately 
equal amounts of energy are radiated by the forward blast wave, moving with 
~ T into the surrounding medium, and by a reverse shock propagating with 
T r — 1 ~ 1 back into the ejecta [39]. The electrons in the forward shock are 
hotter by a factor T than in the reverse shock, producing two synchrotron 
peaks separated by T 2 , one peak being initially in the optical (reverse) and 
the other in the 7/X band (forward) [41,42]. Detailed calculations and pre- 
dictions of the time evolution of such a forward and reverse shock afterglow 
model ( [43]) preceded the observations of the first afterglow GRB970228 ( 
[13,88]), which was detected in 7-rays, X-rays and several optical bands, and 
was followed up for a number of months. 

The simplest spherical afterglow model concentrates on the forward blast 
wave only. For this, the flux at a given frequency and the synchrotron peak 
frequency decay at a rate [43,49] 

Fu0 , t i3-^-m/a+2 9 ) ) (1) 

where g is the exponent of T oc r~ 9 and (3 is the photon spectral energy slope. 
The decay rate of the forward shock F v in equ.(l) is typically slower than 
that of the reverse shock [43], and the reason why the "simplest" model was 



stripped down to its forward shock component only is that, for the first two 
years 1997-1998, afterglows were followed in more detail only after the several 
hours needed by Beppo-SAX to acquire accurate positions, by which time 
both reverse external shock and internal shock components are expected to 
have become unobservable. This simple standard model has been remarkably 
successful at explaining the gross features and light curves of GRB 970228, 
GRB 970508 (after 2 days; for early rise, see §IV) e.g. [96,85,91,73]. 

This simplest afterglow model has a three-segment power law spectrum 
with two breaks. At low frequencies there is a steeply rising synchrotron self- 
absorbed spectrum up to a self- absorption break v a , followed by a +1/3 energy 
index spectrum up to the synchrotron break v m corresponding to the minimum 
energy 7 m of the power-law accelerated electrons, and then a — (p — l)/2 
energy spectrum above this break, for electrons in the adiabatic regime (where 
is the electron energy distribution above r y m ). A fourth segment and a 
third break is expected at energies where the electron cooling time becomes 
short compared to the expansion time, with a spectral slope — p/2 above that. 
With this third "cooling" break z/^, first calculated in [47] and more explicitly 
detailed in [80], one has what has come to be called the simple "standard" 
model of GRB afterglows. One of the predictions of this model [43] is that 
the relation between the temporal decay index a, for g = 3/2 in T oc r~ 9 , is 
related to the photon spectral energy index j3 through F v oc t a v^ , with at = 
(3/2)/3 . This relationship appears to be valid in many (although not all) cases, 
especially after the first few days, and is compatible with an electron spectral 
index p ~ 2.2 — 2.5 which is typical of shock acceleration, e.g. [91,80,95], etc. 
As the remnant expands the photon spectrum moves to lower frequencies, and 
the flux in a given band decays as a power law in time, whose index can change 
as breaks move through it. For the simple standard model, snapshot overall 
spectra have been deduced by extrapolating spectra at different wavebands 
and times using assumed simple time dependences [92,95]. These can be 
used to derive rough fits for the different physical parameters of the burst 
and environment, e.g. the total energy E, the magnetic and electron-proton 
coupling parameters and e e and the external density n . 

IV "POST-STANDARD" AFTERGLOW MODELS 

The most obvious departure from the simplest standard model occurs if the 
external medium is inhomogeneous: for instance, for n ex t oc r~ d , the energy 
conservation condition is Y 2 r 3 ~ d ~ constant, which changes significantly the 
temporal decay rates [47]. Such a power law dependence is expected if the 
external medium is a wind, say from an evolved progenitor star as implied in 
the hypernova scenario (such winds are generally used to fit supernova remnant 
models). Another departure from a simple impulsive injection approximation 
is obtained if the mass and energy injected during the burst duration t w (say 



tens of seconds) obeys M(> T) oc T~ s , E(> T) oc T 1_s , i.e. more energy 
emitted with lower Lorentz factors at later times (but still shorter than the 
gamma-ray pulse duration). This would drastically change the temporal decay 
rate and extend the afterglow lifetime in the relativistic regime, providing a 
late "energy refreshment" to the blast wave on time scales comparable to the 
afterglow time scale [72] . These two cases lead to a decay rate 



Expressions for the temporal decay index a(f3, s, d) in F v oc t a are given by 
[47,72], which now depend also on s and/or d. The result is that the decay 
can be flatter (or steeper, depending on s and d) than the simple standard 
a = (3/2)/?. A third non-standard effect, which is entirely natural, occurs 
when the energy and/or the bulk Lorentz factor injected are some function 
of the angle. A simple case is E Q oc Q~\ T Q oc 6~ k within a range of angles; 
this leads to the outflow at different angles shocking at different radii and 
its radiation arriving at the observed at different delayed times, and it has a 
marked effect on the time dependence of the afterglow [47], with a = a([3,j, k) 
flatter or steeper than the standard value, depending on j, k. Thus in general, 
a temporal decay index which is a function of more than one parameter 



is not surprising; what is more remarkable is that, often, the simple relation 
a — (3/2)/3 is sufficient to describe the overall behavior at late times. 

Evidence for departures from the simple standard model is provided by, e.g., 
sharp rises or humps in the light curves followed by a renewed decay, as in 
GRB 970508 ( [64,66]). Detailed time-dependent model fits [62] to the X-ray, 
optical and radio light curves of GRB 970228 and GRB 970508 show that, in 
order to explain the humps, a non-uniform injection or an anisotropic outflow 
is required. These fits indicate that the shock physics may be a function of 
the shock strength (e.g. the electron index p, injection fraction £ and/or e&, e e 
change in time), and also indicate that dust absorption is needed to simul- 
taneously fit the X-ray and optical fluxes. The effects of beaming (outflow 
within a limited range of solid angles) can be significant [60] , but are coupled 
with other effects, and a careful analysis is needed to disentangle them. 

Prompt optical, X-ray and GeV flashes from reverse and forward shocks, 
as well as from internal shocks, have been calculated in theoretical fireball 
shock models for a number of years [41,42,57,43,81], as have been jets (e.g. 
[38,40,42], and in more detail [76,62,60,77]). However, observational evidence 
for these effects were largely lacking, until the detection of a prompt (within 
22 s) optical flash from GRB 990123 with ROTSE by [2], together with X-ray, 
optical and radio follow-ups [34,21,18,1,12,27]. GRB 990123 is so far unique 
not only for its prompt optical detection, but also by the fact that if it were 





F v <x.t a u^ ,with a — a((3,d,s,j,k,- ■ •) , 
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emitting isotropically, based on its redshift z = 1.6 [34,1] its energy would 
be the largest of any GRB so far, 4 x 10 54 ergs. It is, however, also the first 
(tentative) case in which there is evidence for jet-like emission [34,18,12]. An 
additional, uncommon feature is that a radio afterglow appeared after only 
one day, only to disappear the next [21,34]. 

The prompt optical light curve of GRB 990123 decays initially as oc t~ 2 - 5 
to oc t~ L6 [2], much steeper than the typical oc t~ 1A of previous optical af- 
terglows detected after several hours. However, after about 10 minutes its 
decay rate moderates, and appears to join smoothly onto a slower decay rate 
oc t~ 1A measured with large telescopes [21,34,18,12] after hours and days. The 
prompt optical flash peaked at 9-th magnitude after 55 s [2] , and in fact a 9-th 
magnitude prompt flash with a steeper decay rate had been predicted more 
than two years ago [43], from the synchrotron radiation of the reverse shock 
in GRB afterglows at cosmological redshifts (see also [81,41,42]). An origin 
of the optical prompt flash in internal shocks [43,49] cannot be ruled out, but 
is less likely since the optical light curve and the 7-rays appear not to corre- 
late well [82,21]. The subsequent slower decay agrees with predictions for the 
forward external shock [43,82,49]. 

The evidence for a jet is based on an apparent steepening of the light curve 
after about three days [34,18,12]. If real, this steepening is probably due 
to the transition between early relativistic expansion, when the light-cone is 
narrower than the jet opening, and the late expansion, when the light-cone has 
become wider than the jet, leading to a drop in the effective flux [76,34,49,77]. 
A rough estimate leads to a jet opening angle of 3-5 degrees, which would 
reduce the total energy requirements to about 4 x 10 52 ergs. This is about two 
order of magnitude less than the binding energy of a few solar rest masses, 
which, even allowing for substantial inefficiencies, is compatible with currently 
favored scenarios (e.g. [68,37]) based on a stellar collapse or a compact merger. 

V LOCATION AND ENVIRONMENTAL EFFECTS 

The location of the afterglow relative to the host galaxy center can pro- 
vide clues both for the nature of the progenitor and for the external density 
encountered by the fireball. A hypernova model would be expected to occur 
inside a galaxy in a high density environment n > 10 3 — 10 5 cm -3 . Most of 
the detected and well identified afterglows are inside the projected image of 
the host galaxy [6], and some also show evidence for a dense medium at least 
in front of the afterglow ( [52]). 

In NS-NS mergers one would expect a BH plus debris torus system and 
roughly the same total energy as in a hypernova model, but the mean distance 
traveled from birth is of order several Kpc [8] , leading to a burst presumably 
in a less dense environment. The fits of [95] to the observational data on GRB 
970508 and GRB 971214 in fact suggest external densities in the range of n Q = 



0.04-0.4 cm -3 , which would be more typical of a tenuous interstellar medium. 
These could be within the volume of the galaxy, but for NS-NS on average one 
would expect as many GRB inside as outside. This is based on an estimate 
of the mean NS-NS merger time of 10 8 years. BH-NS mergers would occur in 
timescales ~ 10 7 years, and would be expected to give bursts inside the host 
galaxy ( [8]; see however [19]). In at least one "snapshot" standard afterglow 
spectral fit for GRB 980329 [75] the deduced external density is n Q ~ 10 3 
cm -3 . In some of the other detected afterglows there is other evidence for a 
relatively dense gaseous environments, as suggested, e.g. by evidence for dust 
[74] in GRB970508, the absence of an optical afterglow and presence of strong 
soft X-ray absorption [24,51] in GRB 970828, the lack an an optical afterglow 
in the (radio-detected) afterglow ( [86]) of GRB980329, and spectral fits to the 
low energy portion of the X-ray afterglow of several bursts [52]. One important 
caveat is that all afterglows found so far are based on Beppo-SAX positions, 
which is sensitive only to long bursts t b >,20 s [28]. This is significant, since it 
appears likely that NS-NS mergers lead [37] to short bursts with £&<, 10 s. To 
make sure that a population of short GRB afterglows is not being missed will 
probably need to await results from HETE [26] and from the planned Swift 
[84] mission, which is designed to accurately locate 300 GRB/yr. 

log ( E / eV ) 



FIGURE 3. Spec- 
trum (top) of a hypernova 
funnel model for various 
observer times, showing 
(bottom) the total and Fe 
light curves and equivalent 
widths (Wcth, Mcszaros , 
Kallman & Rees 1999), for 
with R = 1.5 x 10 16 cm, 
n = 10 10 cm~ 3 , and Fe 
abundance 10 2 times solar. 
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The environment in which a GRB occurs may also lead to specific spec- 
tral signatures from the external medium imprinted in the continuum, such 
as atomic edges and lines [5,65,46]. These may be used both to diagnose the 




chemical abundances and the ionization state (or local separation from the 
burst), as well as serving as potential alternative redshift indicators. (In ad- 
dition, the outflowing ejecta itself may also contribute blue-shifted edge and 
line features, especially if metal-rich blobs or filaments are entrained in the 
flow from the disrupted progenitor debris [45], which could serve as diagnostic 
for the progenitor composition and outflow Lorentz factor). To distinguish 
between progenitors, an interesting prediction ( [46]; see also [22,9]) is that 
the presence of a measurable Fe K-a X-ray emission line could be a diagnos- 
tic of a hypernova, since in this case one may expect a massive envelope at 
a radius comparable to a light-day where 7t<1, capable of reprocessing the 
X-ray continuum by recombination and fluorescence. Detailed radiative trans- 
fer calculations have been performed to simulate the time-dependent X/UV 
line spectra of massive progenitor (hypernova) remnants [93], see Figure 3. 
Two groups [67,98] have in fact recently reported the possible detection of Fe 
emission lines in GRB 970508 and GRB 970828. 

An interesting case is the apparent coincidence of GRB 980425 with the 
unusual SN Ib/Ic 1998bw [20], which may represent a new class of SN [29,7]. 
If true, this could imply that some or perhaps all GRB could be associated with 
SN Ib/Ic [90], differring only in their viewing angles relative to a very narrow 
jet. Alternatively, the GRB could be (e.g. [97]) a new subclass of GRB with 
lower energy E 1 ~ 10 48 (Qj/4n) erg, only rarely observable, while the great 
majority of the observed GRB would have the energies E 1 ~ 10 54 (f2 ? /47r) 
ergs as inferred from high redshift observations. The difficulties are that it 
would require extreme collimations by factors 10~ 3 — 10~ 4 , and the statistical 
association is so far not significant [32] . However, two more GRB light curves 
may have been affected by an anomalous SNR (see, e.g. the review of [94]). 



VI CONCLUSIONS 



The fireball shock model of gamma-ray bursts has proved quite robust in 
providing a consistent overall interpretation of the major features of these 
objects at various frequencies and over timescales ranging from the short ini- 
tial burst to afterglows extending over many months. Significant progress 
has been made in understanding both the phenomenology and the physics of 
these obejcts, which may the most widely studied type of black holes sources. 
There still remain a number of mysteries, especially concerning the identity 
of their progenitors, the nature of the triggering mechanism, the transport of 
the energy, the time scales involved, and the nature and effects of beaming. 
However, the collective theoretical and observational understanding is vigor- 
ously advancing, and with dedicated new and planned observational missions 
under way, further significant progress may be expected in the near future. 
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